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Abstract 

Whether molecular states indeed exist in nature has been disputed for a long time. Several new 

resonances have been observed in the recent experiments and they seem to be of exotic structures 

and some of them have been proposed to be molecular states. The very recent observation of 

Z 6 (10610) [(10608.4 ±2.0) MeV] and Z fe (10650)[(10653.2 ± 1.5) MeV] encourages the interpretation 

of multi-quark states. In the Beter-Salpeter (BS) approach, we study the possibility if two heavy 
<— j ■ 

mesons can form a molecular state by exchanging light mesons. Our results indicate that two 
q_[ heavy mesons can form an isospin singlet (I = 0) bound state but cannot form an isospin triplet 

| (1 = 1) when the contribution of a— exchange is reasonably small, i.e. as the coupling of a with 

mesons g a takes the value given in previous literatures. Thus we conclude that the newly observed 
__b(10610) should not be a molecular state, but a tetraquark state instead, at most, the fraction of 
the molecular state in the physical resonance Z&( 10610) is tiny. 
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I. INTRODUCTION 



The naive non-relativistic quark model[l| tells us that a meson contains a quark and an 
anti-quark while a baryon is composed of three valence quarks. However, it is also noted 
that the exotic hybrid, glueball and multi-quark states are not excluded by the SU (3) quark 
model. Until now no any meson or baryon has been confirmed as an exotic state even though 
several resonances are considered to be exotic, namely have different component-structure 
from regular hadrons. There is not any principle to forbid their existence, and in fact the 
discussion about the exotic states has never stopped. For example the exotic structure of 
/o(980) and ao(980) is still under intense dispute |2H6j. Since the conventional qq structure 
cannot well fit experimental data of certain modes0, 0] , some authors proposed that they 
may be molecular states made of two color-singlet mesons or tetraquarks composed of a 
color-anti-triplet diquark and a color-triplet anti-diquark j^-Q], . 

Recently, a series of charmonium-like resonances have successively been experimen- 
tally observed, such as X(3872)Q, X(394O)0, 3^(3940) 0, ^(4430) ± [li and several 
bottomoniun-like states are also discovered by the Belle and Babar collaborations, such 
as Z&(10610) and Zf,(10650) 13J. It is noted that there is almost no room in the regular 
representations of 0(3) <g> SUf(3) ® SU S (2) to accommodate those newly observed resonances. 
Especially these particles Z(4430), Z b (10610) and Z b (10650) with non-zero charge cannot be 

. Concretely 
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understood in the conventional QQ structure if they are real resonances 
the mass of Z& is between the mass of T(4S) and T(5S f ), so it should have one b (quark) 
and one b (anti-quark) but it by no means is a bottomonium because of its electric charge. 
Apparently if it is verified to be a resonance it should be a multi-quark state, namely a 
molecular state or tetraquark should be the preferred choice since they are the simplest 
extension beyond the regular QQ structure. 

Guo et al. explored possible KK bound states in the Bethe-Salpeter framework and found 
that the bound state could exist fill. In the same framework the resonance 5* (5725) was 
considered as a BK molecular state 2(J. Whether two heavy mesons can form a molecular 



state has not been thoroughly investigated in those works [19j, |20(. Thus in this work we 
will try to study the question in the Bethe-Salpeter (BS) approach where the relativistic 
corrections are automatically included. Besides the bound states of two pseudoscalars we 
also explore the bound states which contain one or two vectors. For the aim we need to 
deduce corresponding BS equations. Since in this work we only concern the ground states, 
the orbital angular momentum between two constituent mesons is zero (/ = 0). For a 
molecular state whose constituents are one pseudoscalar and one vector, its J PC is 1 + . For 
the molecular states which consist of two vector mesons their J PC may be + , 1 + and 2 + . 
We only deduce the BS equation for the bound state with J p = + because it should be the 
lightest molecular state and more favorable in the nature. As a matter of fact, the Lorentz 
structures of the others are very complicated and will be explored in our later works. 
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When we solve the BS equations for these molecular states, two approximations are 
needed: the ladder approximation and instantaneous approximation. In general the en- 
ergy exchange between the constituents is small compared to A.qqd which characterizes the 
binding energy scale of the constituent mesons, so that the instantaneous approximation is 
applicable. The simple ladder approximation has been employed all along in the history, 



but recently some works suggest that the cross-ladder should be included [2l|, |22|. Since our 
goal is to study the possibility of formation of molecular states instead of making precise 
theoretical predictions on the spectra we will still use the simple ladder approximation to 
gain our qualitative conclusion. 

After this introduction we derive the BS equations for the + and 1 + molecular states. 
Then in section III we present our numerical results of the binding energies along with 
explicitly displaying all input parameters. Section IV is devoted to a brief summary. 



II. THE BETHE-SALPETER FORMALISM 

In this section we will deduce the BS equations for the + and 1 + molecular bound states. 

A. The bound state (0 + ) composed of two pseudoscalar mesons 



In Ref. [19|, |20j the BS equation of a bound state of two pseudoscalars was deduced. 



Since in this work we need to use the corresponding formulas for analyzing the spectra of 



the molecu 
of Ref. E 



ar states containing two pseudoscalars, let us briefly review the main contents 



201 ] . then in the following subsections, we will step forward to deduce the BS 
equations for one-vector-one-pseudoscalar bound states and two-vector bound states . 

The BS wave function for the bound state \V) of two pseudoscalar mesons can be defined 
as following: 

X v {x u x 2 ) = (0|T&(z 1 )&(z 2 )|7>) = e~ wx Xv {x) , (1) 

where 0i(xi) and 02(^2) are the field operators of two mesons, respectively, P denotes 
the total momentum of the bound state, the relative coordinate x and the center of mass 
coordinate X are 

X = 771X1 + 772X2 , x = xx—x 2 , (2) 

where 77, = m;/ (mi + m 2 ) and rrii (i = 1, 2) is the mass of the 7-th constituent meson. The 
equation for the BS wave function can be derived from a four-point Green function, 

S(x 1 ,x 2 ;y 2 ,y 1 ) = (0|T 1 (x 1 )02(x 2 )(0i(z/i)0 2 (y2)) t |O) . (3) 



3 



To obtain the BS equation, we express the above four-point Green function in terms of the 
four-point truncated irreducible kernel K , 



+ / d 4 u 1 d 4 u 2 d 4 v 1 d 4 v 2 S7 ) (x u x 2 ;u 2 , u 1 )K(u 1 ,u 2 ; v 2 , vi)S(v 1 ,v 2 ; y 2 , Vi) , (4) 



S(xi,x 2 ;y 2 ,yi) = S^(xi,x 2 ;y2,yi) 

S (0)l 

where S(p) is related to the forward scattering disconnected four-point amplitude, 

S i0 )(xi,x 2 ] y 2 , yi) = Ai(xi, yi)A 2 (x 2 , y 2 ) , (5) 
and Ai(xi, y^ is the complete propagator of the i-th meson. Here we have 

Aj (p 1 ,m 1 )A 2 (p2,m 2 )x T (p) = J -—-^K^p^x^p') , (6) 



where A a - A 2 - ^^ang^ K i [p, p ) - -i c f g (pi- Pi ) 2 -m* 



2QJ. 



which has been deduced in Ref. 

The relative momenta and the total momentum of the bound state in the equations are 
defined as 

P = V2P1 - V1P2 , P = V2P1 ~ V1P2 > p = Pi + P2 = Pi + P2 ■ ( 7 ) 

Directly solving the BS equation ([6]) is extremely difficult. In general one needs to use 
the so-called instantaneous approximation: Ki(p,p') = ^(p, p') by which the BS equation 
can be reduced to 

£r + ( S^ - (p) = U w? k ^( P , P ')xm)f( P - p') 2 , (8) 



where Ei = yp 2 + mf, E = P°, and the equal-time wave function is defined as X P (p) = 
J dp x P {p) ■ Since the constituent meson is not a point particle a form factor at each in- 
teraction vertex among hadrons must be introduced to reflect the finite-size effects of these 
hadrons. The form factor is assumed to be in the following form: 

where A is a cutoff parameter. For exchange of a light vector between the mesons, the kernel 
is 

v- i i\ ■ i (p + p') 2 + A VlV2 E" + (p 2 - p' 2 ) 2 /M 2 

K iv (p,p) = ic l9pvp g pvp (p _ p /)2 + M 2 • ( 10 ) 

The exchanged mesons between two pseudoscalars are vector mesons, obviously we only 
need to keep the lightest vector mesons p and u for taking the dominant contributions into 
account |l9|, |20]. When the bound state is an isospin-scalar c/ = 3 for p (n) and cj = 1 for 
uj (a). When the bound state is an isospin- vector ci = —1 for p (ir) and cj = 1 for u (a). 
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B. The bound state(l + ) composed of a pseudoscalar and a vector 

We can define the BS wave function for the bound state |V) composed of one pseudoscalar 
and a vector mesons as following: 

X\(xi,x 2 ) = (0|T&(xi)#(x a )|V) = e~* P V», (11) 

where <pi(xi) and (f> 2 (x 2 ) are respectively the field operators of the two mesons. The equation 
for the BS wave function should be derived from a four-point Green function, 

S^(x 1 ,x 2 ; y 2 , Vl ) = (0\T^( Xl )Mx2)^i(yi)My2)T\0) • (12) 

To obtain the corresponding BS equation, we write the above four-point Green function in 
terms of the four-point truncated irreducible kernel K , 

S^ixu x 2 ; y 2 , Vi) = S^(x u x 2 ; y 2 , y x ) 

+ J d^uid^u^vid^ Sffi(xi, x 2 \ u 2 , Ui)K a p{u l} u 2 ; v 2 , v^S^iv^ v 2 ; y 2 , y x ) , (13) 

where S(p) is related to the forward-scattering disconnected four-point amplitude, 

Sj^ixuXz; V2, Vi) = Ai(xi, y 1 )A 2 ia (x 2 , y 2 ) , (14) 

and Aj(xj, yi) is the full propagator of the i-th particle, 

^\ Pl , mi )A^ a (p 2 ,m 2 ) X \(p) = J -0^K*^p')x" v (P) ■ ( 15 ) 

Here A 1 = p ^ m 2 and A 2/ia = ^1^2 {^r 1 - g m ) are the propagators of pseudoscalar and 
vector mesons. 

From Eq. (fl~5l) one can obtain 



iA - m i^- mi f M = / $k«MAV)<!^ - <n • («> 

We write x fl v {p) = X v (p) e/ \ then multiply an e* on both sides. Summing over the 
polarizations we deduce a new equation 

T-2 2^2 T.XM = I ^K a ,(p,p') Xv (p')(^4 - gn^K- ~ 4) ■ (17) 

[Pi — mi)(p 2 — m 2 ) J {2np v m 2 p z p 

With the Feynman diagrams depicted in Fig.l and the effective interactions (13. fl^ we even- 
tually obtain 

_ ^f-g XlT 

Ka/3 = g PVP g vvv — 2 jJ72~(Pl + Pl)\[-g a /3{P2 + P 2 )o + P 2a g*i3 + P^&ra]- 

q — My 



Ma ' f/3 

Defining K 2 (p,p') = K a ^(p,p')(^r- — and employing the instantaneous 

approximation we obtain the BS equation similar to Eq. ([6]) and ([8]) but with a different 
kernel, 



X V (P) = \ J K 2V(S)(P,P')X V (P')F{P ~ P') 2 



(E 1 + E 2 )/E 1 E 2 /Xv ^ J 2 J (2tt 



j pi 



lh 



P2 




(a) (b) 

FIG. 1: (a) A bound state composed of two pseudoscalars (b) a bound state composed of a pseu- 
doscalar and a vector. In the two Feynman diagrams p and uj are exchanged. 



f Pl q' A i Pi v pi a , n pi 




P2 y „ P2 



(a) (b) 
FIG. 2: A bound state composed of two vectors, (a) p and oj are exchanged (b) ir is exchanged 



C. The bound state (0 + ) composed of two vector mesons 

The quantum number J p of the bound state composed of two vectors can be + , 1 + and 
2 + . As aforementioned, since + is more favorable in nature and its Lorentz structure is 
relatively simpler than the two others, in this work we only study the + bound states and 
define the corresponding BS wave function \S) as following: 

X s {xi,x 2 ) = (0|T^(xi)<^(x 2 )|5) = e- wx Xs (x). (19) 

The equation for the BS wave function can be derived from a four-point Green function, 

S( Xl ,x 2 ;y 2 , yi ) = (OlT^^)^^)^^)^^))^). (20) 
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In analog to the procedures used in last subsection, we obtain 



AlilaiPl,™!)^ 1 ^' (P2,m 2 )X s (p) = I 77r^ K a (PiP')X s (p') 



where A 



Ilia „2_ m 2 



(27T) 4 ' 



(21) 



, i (Pi*P^-g ) an d A 



From Eq. (fToT) one can obtain 
-1 



X%(p) 



(27T) 



(p,^)xs(?/)(^ - <n( £1 ^ - (22) 

r«2 r«2 



19 



20 
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With the Feynman diagrams depicted in Fig. 2 and the effective interaction 
we can obtain 

Kv* = 9vw9' vvv — 9 , — ^l-g a 'Hpi+p[)\+p[ a '9\+p^g\]l-g a t3(P2+p' 2 )*+P2 a 9*p+P2pgo 
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'VVV q 2 _ M 2 



and 



K Pa = gvpvg' VPV £a Pr7P q*Plp£af3aipiq a P2 ■ 



Defining K 3V(PtS) (p,p') = K~ {m Jp,jf){^ - <r0(*^ - <w) we derive the BS 
equation which is similar to Eqs. fl6]) and © but possesses a different kernel, 



E 2 -(E l + E 2 ) 2 ^ , , i 



{E 1 + E 2 )/E 1 E, 



■XAP) 



d 3 p' 



2 J (2tt) 



#3V(p,s)(p, P / )x 5 (p / ) ir (P - P 



f\2 



(23) 



III. NUMERICAL RESULTS 



Solving the BS equations Eq. OS), PSD an d we obtain the eigenvalues for each 

state. To determine whether such a state composed of two heavy mesons is a bound state, a 
criterion must be set. That is: if the eigenvalue of the bound state obtained by solving the 
BS equation is negative, namely the total mass of the system is lower than the sum of the 
masses of the two constituents, such a system is considered as a bound state, i.e. a hadronic 
molecule which may exist in the nature. 

Now let us solve the BS equations. Since the function Xv(p), Xv(p), Xs(p) only depends 
on the norm of the three-momentum we may first integrate over the azimuthal angle of the 
function in (B]h (ITS]) or (125]) 

i/g*(p,p W p-pr, 

to obtain a new form f7(|p|, |p'|), then the BS equation turns into a one-dimension integral 
equation 
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Letting |p|(|p'|) take ordered n discrete values and the gap between two adjacent values 
be Ap, then x(|p|) can be arranged as a column matrix and the coefficients constitute an 
n x n matrix M. The scalar equation of M would determine the n eigenvalues E n and the 
smallest one should correspond to the binding energy of the ground state. In our calculation 
we take n = 151 and the first discrete value of |p| = 0.001 GeV and the largest |p| is 2 GeV. 

In our calculation we need to determine the values of the parameters A, g PVP , g vvv , g PSP 
and g vsv - By fitting the data of the decay processes where the couplings KKp andBBp are 
concerned, the authors of Ref. 20J obtained the cutoff parameter A which can vary from lGeV 
to 4GeV. Their strategy was to fit the measured decay widths in terms of the theoretically 
derived formulas with A, thus they obtained the rather wide range for A. Indeed such a wide 
range reduces our theoretical prediction power. Instead, in Ref. 25| the authors suggested 
a relation: A = m + a Aqcd where m is the mass of the exchanged meson, a is a number 
of 0(1) and Aqcd = 220 MeV i.e. A ~ lGeV for exchanging p or uj. Since there are 
uncertainties on A we will choose A = 1.5 GeV, A = 2.5 GeV and A = 3.5 GeV to calculate 
the eigenvalues and compare the results. The coupling constants gsB P = gi<K P = 3 are fixed 
by fitting data 23j. Considering the flavor SU (3) symmetry and heavy quark spin symmetry 

we believe that g B B P (u) = gDD P (w) = Sbb* p (w) = 9dd* p {w) = gB*B* P (oj) = gD*D* P (w) = g should 
be a not-bad approximation. Instead of setting g = 3, we also let g take various values from 
2 to 4 to investigate dependence of the eigenvalues on those coupling constants. 

The masses of the concerned constituent mesons are directly taken from the databook 



26j as: m D = 1.865 GeV, m D » = 2.007 GeV, m B = 5.279 GeV and m B * = 5.325 GeV 



A. The exchanged mesons are p and u) 

It is noted that in those tables there are many places symbolized by "— " which means such 
bound states cannot exist because either the BS equation with the corresponding parameters 
has no solution at all, or the derived eigen-energy is not negative, namely the mass of the 
supposed-to-be bound state is larger than the sum of the two constituent mesons. 

We calculate the eigenvalues of these bound states of DD(0 + ), DD*(1 + ), D*D*(0 + ), 
BB(0 + ), BB*(1 + ), B*B*(0 + ) respectively. Apparently when the parameters A and 
9pvp(.9vvv) are reasonable, the corresponding BS equation is solvable. For the bound states 
of isospin / = we can obtain the ground eigenvalues (see Tab. HI). It implies that the two 
heavy mesons can form an isospin-0 molecular state by exchanging light mesons. However 
for / = 1 we cannot obtain eigenvalues which are below the threshold with the same param- 
eters. The reason is that the contributions of exchanging p and u nearly cancel each other 
in the kernel K. 
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TABLE I: The eigenvalues of the ground molecular states (I = 0) when only p and u are exchanged 
between the constituents 





DD 


DD* 


D*D* 


BB 


BB* 


B*B* 


A — 


1.0 




3.877 


3.406 


10.50 


10.54 


9.904 


A = 


2.5 a 


3.708 


3.844 


2.699 


10.38 


10.43 


9.245 


A = 


3.5 a 


3.680 


3.815 


2.350 


10.33 


10.37 


8.950 


A = 


1.5 6 


3.618 


3.756 


3.892 


10.26 


10.31 


8.899 


A = 


2.5 6 


3.379 


3.529 


3.845 


9.955 


10.00 


7.771 


A = 


3.5 6 


3.254 


3.414 


3.841 


9.800 


9.857 


7.283 


A = 


1.5 C 


3.330 


3.480 


3.637 


9.909 


9.956 


7.798 


A = 


2.5 C 


2.860 


3.067 


3.878 


9.392 


9.446 


6.356 


A = 


3.5 C 


2.656 


2.886 


3.935 


9.162 


9.219 


5.802 



a the coupling constants g PVP{vvv) are set to be 2 
b the coupling constants g PVP(vvv) are set to be 3 
c the coupling constants g PVP{vvv) are set to be 4 



TABLE II: The eigenvalues of the ground molecular states when the contribution of 7r are included 
with g D »j=)»^ = gB*B*Tx = 18 and the coupling constants 9 PVP(VVV) are set to be 3. 





D*D*(I = 0) 


B*B* (I = 0) 


D*D*(l=l) 


B*B*(I = 1) 


A = 


1.5 


3.588 


8.489 






A = 


2.5 


3.606 


7.270 






A = 


3.5 


3.640 


6.760 







B. As exchange of tt between two vectors is taken into account 

For the bound states composed of two vectors such as B*B* and D*D*, the contribution 
of pion-exchange should be included. Including tt— exchange in the kernel, we repeat our 



derivation and obtain a new BS equation. Using the data of D* + — > D°ir + 26( and the formula 



given in Ref. [27J we fix go'Dw = 18; an d furthermore under the heavy quark symmetry we 
may set gB*B*n = 9d*d*tt = 9d*d-k- With this coupling constant we calculate the eigenvalues 
of the bound states D*D* and B*B* and the results are listed in TablLTI Comparing the 
values in TablLTI with that in Tab]TJ we find that the masses of the molecular states D*D* 
and B*B*(I = 0) slightly decrease by less than 10% as the pion-exchange between the two 
vectors is taken into account. It shows that the pion-exchange contributes a weak attractive 
effect, but it is not as important as that from the vector exchange. Moreover, for the 1=1 
bound state the pion-exchange with Cj = — 1 is even less important so there is still no / = 1 
bound state as indicated by our numerical results. 
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TABLE III: The masses of the ground molecular states when one set A = m + oAqcd with 
9bb p = 9bBuj = 3 and g D *D* n = 9b*b*-k 





BB 


BB* 




BB 


BB* 






(1 = 0) 


(1 = 0) 


(1 = 0) 


(1=1) 


(1=1) 


(i = i) 


a = 1 


10.49 


10.53 


9.939 








a = 2 


10.39 


10.44 


9.364 









a the contribution of exchanged 7r has been included. 

TABLE IV: The masses of the ground molecular states and the coupling constants 9 pv p(vw) 
set to be 3 when the contribution of a are included. 





DD 


DD* 


D*D* a 


BB 


BB* 




A = 2.5, g a = 0.76(1 = 0) 


3.386 


3.538 


3.839 


9.967 


9.990 


7.750 


A = 2.5, g a = 0.76 (I = 1) 














A = 2.5, ^ = 3(1 = 0) 


3.505 


3.382 


3.750 


10.15 


9.801 


7.130 


A = 2.5, g a = 3(I = l) 






3.760 




10.56 


10.22 



a the contribution of exchanging ir is not included. 

As aforementioned one also can set A = m + oAqcd which would determine values of 
A lower than that we have used and the relation implies different A values for different 
exchanged mesons. We calculate the masses of the _£?(*)_£?(*) system for a = 1 and a = 2 
respectively and the numerical results are presented in Tab JIIII The results in Tab JIIII are 
qualitatively consistent with the conclusions we have made. 



C. As exchange of a is taken into account 

Moreover, two pseudoscalars and two vectors can couple to a scalar, so that one should in 
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elude the contribution of a— exchange in the kernel. In terms of the effective Lagrangian 
we deduce the corresponding kernel (See appendix). With the contribution of a the eigen- 
values of the bound state DD(0+), DD*(1+), D*D*(0 + ), BB(0 + ), BB*(1+), B*B*(0+) are 



re-calculated. The parameter g a is determined by the 
ory (HHChPT) as g a = ^ « 0.76 with g T = 3.73 1 



leavy hadron chiral perturbation the- 
. Using this value of g a we repeat our 



computations and the resultant masses of the ground molecular states are listed in Tab JIVI 
Comparing TabJIwith Tab il VI we find the contribution of a is very small for the concerned 
case. 



1 We thank Dr. Valery E. Lyubovitskij for explaining how to determine the value of g^. 
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The above numerical results show that one cannot obtain an / = 1 bound state which 
is composed of two heavy mesons. The reason is the contributions from exchanges of p 
and u cancel each other for the quantum number 1 = 1, and the small coupling constant 
g a determines that including the a contribution cannot make a substantial change to the 
kernel. Let us take another angle to look at this issue. With the understanding, now let us 
deliberately enhance the contribution of a by enlarging its coupling with pseudoscalar and 
vector mesons. Namely, we investigate to what value the coupling constant g a reaches, the 
two heavy mesons can form a molecular state. We notice that as the coupling constant g a = 3 
the eigenvalues for BB*(I = 1, J p = 1 + ), B*B* (I = 1, J p = 0+) and D*D*(I = 1, J p = 0+) 
meet the criterion for forming stable molecular states and the results are shown in Tab JIVI 
Moreover, as we set the coupling constant g a = 4 and A = 2.5 we can also obtain the mass of 
3.864 GeV for DD*(I = 1, J p = 1 + ). It is worth emphasizing that when we set A = 2.5 GeV 
and g a = 3 the mass of the bound state BB*(I = 1, J p = 1 + ) is 10.56 GeV which is close to 
the mass of the newly observed narrow-structure resonance Zfc(10610) [(10608. 4±2.0) MeV]. 



IV. A BRIEF SUMMARY 

In this work we study the possibility of two heavy mesons forming a hadronic molecule. 
We employ the BS framework to approach this goal because it may include the relativistic 



corrections automatically. In Ref. 20] the BS equation for the bound state of two pseu- 
doscalar mesons was deduced whereas in this work we derive the BS equations of the bound 
states composed of one vector and one pseudoscalar with J p = 1 + and the bound states 
composed of two vectors with J p = + . Our numerical results indicate that when the pa- 
rameters are within reasonable ranges two mesons indeed can form an / = molecular state 
but cannot form an / = 1 molecular state because the contribution of a— meson is small, 
thus we conclude that the newly observed Zft(10610) prefers a tetraquark structure or at 
most has a tiny fraction of the molecular state of BB*. By contrast, if the coupling constant 
g a is enhanced by several times an I = 1 molecular state might be formed. For example 
when we set A = 2.5 GeV and g a = 3 the mass of the 1 = 1 bound state BB*(1 + ) is 10.56 
GeV which is close to the mass of the newly observed narrow-structure resonance Zb(10610). 

Therefore our qualitative conclusion is that either the observed ^(10610) is not a molec- 
ular state or the coupling of a— meson with the pseudoscalar and vector mesons is several 
times larger than that given in the HHChPT. 

Since the parameters are fixed in various experiments and can span a relatively large 
range we cannot expect all the numerical results which depend on the parameters, to be 
very accurate. The goal of this work is to study the possibility of two heavy mesons can 
form a molecular state. Our results, even not accurate, have obvious qualitative significance. 

In a recent work, Ali et al. [30| propose the tetraquark interpretation of the charged 
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bottomonium-like states Z b (10610) and Z b (10650), and find that with this ansatz their 
estimates on the production and decay fit the data well, thus they claim that the tetraquark 
interpretation is supported by the experimental measurments. This is consistent with our 
conclusion that the tetraquark structure of (10610) is preferred. Definitely, further the- 
oretical and experimental works are badly needed for gaining better understanding of the 
exotic structures of the heavy mesons. 
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C MMp = ig M M P [M^(d^M) - (^M)t^M] • ^ 
Cmm. = iguMuWid^M) - (<9 M M)M] • ^ 

£ p m*m* = ig P M*M*l^ ■ i-M*»l>(t, - t,)M*J + M* v ^d v M*J - PM^-tM^], 
m M * 



C m *m*c = 2m M *g a M* • M*V 

C-mmv = -2m M g*M ■ M^a (Al) 
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where gi, g[, g2, g' 2 , #3 and g' 3 represent the coupling constants in the vertexes and 
q 2 = (p 2 + p ' 2 - 2p • p'). 
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